Charge-induced spin polarization in non-magnetic organic molecule Alq 3 
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Electrical injection in organic semiconductors is a key prerequisite for the realization of organic 
spintronics. Using density-functional theory calculations we report the effect of electron transfer 
into the organic molecule Alq3. Our first-principles simulations show that electron injection sponta- 
neously spin-polarizes non-magnetic Alq3 with a magnetic moment linearly increasing with induced 
charge. An asymmetry of the Al-N bond lengths leads to an asymmetric distribution of injected 
charge over the molecule. The spin-polarization arises from a filling of dominantly the nitrogen p z 
orbitals in the molecule's LUMO together with ferromagnetic coupling of the spins on the quinoline 
rings. 

PACS numbers: 87.15.-v,85.75.-d,70 



Organic semiconductors (OSCs) are contemporary be- 
ing considered as prospective materials for spintronics 
applications [HQ- As OSCs primarily consist of 7r- 
conjugatcd molecules composed of atoms with low atomic 
numbers, the spin-orbit and nuclear hyperfine interac- 
tions are weak. OSCs therefore exhibit typically very 
long spin relaxation times a property favorable for 
spin transport. Moreover, owing to the availability of 
a large variety of organic semiconductors (OSCs) and 
the possibility to chemically tailor specific functionalities, 
there is a huge potential for the development of versatile, 
low-cost organic spintronic devices [3, HI]. 

A key pre-requisite for the realization of organic spin- 
tronics is the realization of electrical injection in the 
organic layer. Recently, the injection of spin-polarized 
carriers from a ferromagnetic contact into an OSC has 
been demonstrated @, [J- Notwithstanding this affirma- 
tive observation, not all OSC materials are comparably 
well suited for organic spintronics. In particular, the 
7r-conjugated OSC tris(8-hydroxy-quinoline) aluminum 
(Alq3) has been discovered to perform very good in or- 
ganic spin-valves, providing one of the best giant magne- 
toresistance (GMR) values Alq3 is a relatively well- 
known material that is widely used as component in or- 
ganic light-emitting diodes (OLEDs) d®. 

However, in 

spite of this, it is not understood why precisely Alq 3 gives 
rise to extremely long spin- relaxation times [3| , to a high 
GMR 0, G3, and significant tunnel magnctorcsistance 
(TMR) values [n|. 

Here we employ density functional theory (DFT) calcu- 
lations to explore the influence of electron injection into 
Alq3. Our first-principles calculations reveal that elec- 
tron charging causes, unexpectedly, a spin-polarization 
of the organic molecule, which is non-magnetic in its neu- 
tral state. The origin of the appearing spin-polarization 
is analyzed in detail and shown to be connected to the 
small asymmetry of the 7r-conjugated rings, i.e., a feature 
that is closely related to the particular molecular geome- 
try and molecular orbitals of Alq 3 . The here-reported in- 
fluence of electron injection in Alq3 provides an essential 



ingredient to a microscopic understanding of the complex 
spin-polarized transport in Alq3. 

Since Alq 3 is widely used as a material in organic 
electro-luminescent devices [H, a number of exper- 
imental as well as theoretical studies have been per- 
formed for this material. Using DFT calculations, photo- 
emission, and near-edge x-ray absorption fine structure 
(NEXAFS), Andreoni and co-workers [l2[ have charac- 
terized the structural and electronic properties of Alq 3 in 
both neutral and charged states, which provided a pic- 
ture of the orbital structure of Alq 3 . Johansson et al. 



13| used a combination of x-ray and ultraviolet photo- 



emission spectroscopy (UPS) with DFT calculations to 
discuss the interaction of Alq 3 with Li and K. Recently, 
the geometrical and electronic structures of Alq 3 inter- 
acting with Mg and Al have also been investigated using 
DFT calculations 14, 15 1. Baik et al. [l6| tried to make 
this material magnetic by Co doping. Zhan et al. [i~7j 
studied the interaction of Alq 3 with an Fe surface. Us- 
ing DFT, semi-empirical, and ah initio molecular orbital 
theories, Zhang et al. [18[ studied the effect of electrical 
charging on the electronic structure of Alq 3 . However, 
these studies have not considered any geometric struc- 
tural change in the molecule due to charge injection, nor 
any magnetic effects. 

We have performed first-principles DFT calculations 
using the full-potential Vienna ab initio simulation pack- 
age (VASP) [19( which uses pseudopotentials and the 
projector augmented wave approach. We used the 
Perdew-Burke-Ernzerhof generalized gradient approxi- 
mation (GGA) for treating the exchange-correlation po- 
tential. A sufficiently high energy cutoff (550 eV) was 
used in each calculation to obtain accurate results. In 
our simulations the forces on each of the atoms were cal- 
culated using the Hcllmann-Feynman theorem, and were 
subsequently used to perform a conjugate gradient struc- 
tural relaxation. The structural optimizations were con- 
tinued until the forces on the atoms converged to less 
than 1 meV/A. This optimization has been completely 
carried through for each Alq 3 molecule with an additional 
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FIG. 1: (Color online) Geometry and charge-induced mag- 
netization density of meridianal Alq;?. Al, O, N, C, and 
H atoms are shown as magenta, red, blue, yellow, and 
turquoise spheres, respectively. The magnetization density of 
the charged Alq3 molecule with one extra electron is shown 
by the green hyper-surface. 



TABLE I: Ground-state Al-N and Al-0 bond lengths (in 
A), for neutral meridianal Alq3 from different theoretical cal- 
culatio ns [l3l . l2(il . |22| , including this work, and experimental 
results [2y] (averaged over several available data sets). 





HF/ 


B3LYP 


B3LYP 


Exp. 


GGA 


bond 


3-21+G** 


[20J 6-31G(d)[22] DNP[13J 


ave.[20] 


our cal. 


Al-Ni 


2.063 


2.08 


2.019 


2.044 


2.063 


Al-N 2 


2.033 


2.06 


2.019 


2.025 


2.050 


AI-N3 


2.117 


2.13 


2.051 


2.079 


2.103 


Al-Oi 


1.826 


1.86 


1.839 


1.847 


1.865 


AI-O2 


1.866 


1.89 


1.868 


1.867 


1.888 


AI-O3 


1.856 


1.88 


1.864 


1.856 


1.886 



TABLE II: Ground state Al-N and Al-0 bond lengths for 
meridianal Alq3 upon injection with 0.5 and 1.0 electron, re- 
spectively. 

bond 0.5 e (in A) 1.0 e (in A) 

Al-Ni 21576 2T073 

AI-N2 2.042 2.031 

AI-N3 2.065 2.048 

Al-Oi 1.877 1.887 

AI-O2 1.899 1.905 

AI-O3 1.889 1.896 



charge. There are two isomers of Alq3, namely facial 
and meridianal with C3 and Ci point group symmetries, 
respecti vely . Of these two, the meridianal isomer is re- 
ported 12 , [2^ to possess the lowest formation energy. 
This we have also confirmed in our calculations. In the 
following we consider therefore only meridianal Alq 3 . 

The molecular structure of the meridianal isomer is 
shown in Fig. Q] The Al atom is surrounded by three 
nearly planar nitrogen atoms as well as by three nearly 
planar oxygens, being the link to the quinolines that con- 
sists of a pyridine and a benzene ring. In our study we 
first adopted the experimental structure (see, e.g., (2p| |) 
and performed a full geometric relaxation for the neutral 
ground state. Subsequently, we used this DFT-computed 
ground state geometry and added extra charges (indi- 
cated later on by q) on the molecule to simulate the 
charge injection and repeated the full atomic relaxation 
procedure, assuming an equivalent neutralizing back- 
ground in the simulation cell. Note that the added charge 
is not set to have any fixed spin. Since in the charged 
molecule there will always be a net electric dipole mo- 
ment, the dipole correction has also been included in our 
calculation. We have also tried positive charging, i.e., ex- 
tracting electron charge from the molecule. However, we 
find that in this case the molecule become very unstable 
and we are not able to converge to a stable molecular 
geometry. This finding is in fact consistent with an ex- 
perimental investigation that shows a very high positive 
charge injection energy as compared to a low energy for 
negative charge injection 23 1. 

To start with, we compare our geometry optimized 



ground state structure of Alq 3 with available theoretical 
and experimental data, which are summarized in Table 
HJ Previous studies showed that structural parameters 
of meridianal Alq 3 in good agreement with experiment 
could be obtained using Hartree-Fock calculations with 
the 3-21+G** basis set [20(. A somewhat less good agree- 
ment was achieved with hybrid-functional B3LYP calcu- 
lations 13 , [22f ■ It can be observed from Table U that 



our GGA calculations yield structural results which are 
in very good agreement with previous Hartree-Fock [ItJ 
and recent B3LYP calculations 22 1 . 

One observes from the structural data that the Al is 
bonded more closely to the oxygens than the nitrogens, 
and that it has different bond lengths to the six ligands. 
The N3 atom has the largest separation from Al whereas 
Ni and N2 have similar bond lengths with Al. The three 
Al— O bond lengths are more equal. Given the different 
Al— N bond lengths in an otherwise symmetric arrange- 
ment, an influence of charge doping on the molecule is 
not surprising. It is therefore interesting to examine the 
behavior of the molecule upon electron charging. 

Table |H] shows the structural change of the molecule 
due to electron charging. Compared to the charge- 
neutral structure, the Al— N2, Al— N3 bond lengths are 
reduced, Al— Ni becomes slightly longer. The Al— Oi 
bond length also becomes longer, but the Al— O2 and 
Al— O3 distances are not changed as much. The change 
is gradual with the amount of charging of the molecule. 
The influence of the charge is most prominent for the 
N3 atom having initially the maximum bond length. We 
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FIG. 2: (Color online) Magnetic spin moment of Alq3 as a 
function of injected electron charge. In the inset, the pro- 
jected moments on the thee N atoms are shown with the added 
electron charge. 

further observe that the asymmetry of the three Al— O 
bond lengths is reduced with charging. 

Our calculations reveal that upon electron charging 
the non-magnetic Alq3 molecule spontaneously spin- 
polarizes. Fig.[2]shows the development of the molecule's 
spin moment due to the added charge. The magnetic 
moment of the molecule increases approximately linearly 
with the extra charge. This is a surprising result, be- 
cause the extra charge could be equally distributed over 
the two spin channels of a molecular orbital, leading to 
a non-magnetic charged molecule, but this does not hap- 
pen. A key to understanding why this is the case is 
the asymmetric bonding in the neutral molecule. The 
added electron charge docs not distribute equally over 
the three nitrogens, but it goes dominantly to the N3 
atom. In Fig. [T] the magnetization density corresponding 
to q = 1.0 e is shown by the hyper-surface. It is clearly 
observed that the magnetization density resides promi- 
nently on the three N atoms. For this relatively high q 
value, the carbon atoms of the pyridine ring display also 
some magnetization density, corresponding to magnetic 
moments of 0.08 fiB on the average. This magnetization 
density resides mainly on one of the carbon atoms near- 
est to nitrogen, as well as on the 3 rd nearest neighbor 
carbon atom. Neither O nor Al atoms carry a magnetic 
moment. 

From the above we deduce that the magnetic moment 
arising due to electron charging is related mostly to the 
N atoms. To understand this further, we have plotted 
the local magnetic moments on the three nitrogens as a 
function of extra electron charge on the molecule. This 
is shown in the inset of Fig. [2j The magnetic moments 
on the Ni and N2 atoms are varying in a similar man- 
ner whereas the curve for N3 is separated from that of Ni 
and N2 . Over the full range of charging, the N3 atom has 
acquired a larger magnetic moment than the other two 
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FIG. 3: (Color online) Atom projected density of states 
(DOS) of the Alq3 molecule shown in (a) for the non- 
magnetic, charge neutral, and in (b) for the molecule with 
1.0 electron doping. In (b), the spin-resolved partial DOS is 
shown. The arrow indicates the shift of the LUMO to the 
Fermi energy (E F ). Note that the N states of the LUMO 
become spin-polarized upon charging. 



N atoms. This behavior corroborates with the structural 
change upon charging. As discussed earlier, the decrease 
of the Al— N3 bond length is more significant than that 
of Al— Ni and Al— N2. This larger bond length reduction 
appears to be related to the steeper increase in the mag- 
netic moment of N3. The Al— O bond lengths are not 
modified appreciably with extra charge and consistently 
no spin-polarization on the oxygens is observed. 

The atom-projected density of states (DOS) are shown 
in Fig. [3] for both neutral and a molecule charged with 
q = 1.0 e. The highest Al p states occur at 2 eV bind- 
ing energy. The Al orbitals are hence saturated with 
covalent bonds to the oxygens and nitrogens. The Al 
3p shell is in fact over-saturated in the six-fold coordi- 
nation, which explains why one nitrogen, here N3, has 
a longer distance to Al. The quinoline ring is however 
structurally very stable, giving not much spacial freedom 
for a longer Al— N distance. As can be seen from Fig. 
[3ji for the neutral molecule, the HOMO level has mainly 
O-p character whereas the LUMO consists of mainly N-p 
states with a small admixture of O and C states. Al has a 
negligible contribution to both HOMO and LUMO. The 
calculated HOMO-LUMO gap is 1.87 eV. Upon charging, 
the N-p states in the LUMO become spin-polarized and 
the HOMO-LUMO gaps for the spin-up and spin-down 
channels increase to 2.0 and 2.2 eV, respectively. The 
HOMO states remain practically non-polarized. 

The origin of the unexpected magnetization in the 
electron-injected state can be understood from the fol- 
lowing observations. Due to the molecule's asymmetry 
the unoccupied N p-dominated states at 1.8 eV occur at 
slightly different energies, with the N3 related states oc- 
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FIG. 4: (Color online) Atom-resolved p partial density of 
states (DOS) of the lowest unoccupied molecular orbitals of 
Alq3, shown separately for the three quinoline groups (labeled 
in Fig.[l|. The p states of the quinoline-3 farthest away from 
Al have the lowest energy. 



curring at a somewhat lower energy. Fig. |3] shows the un- 
occupied p partial DOS of the atoms in the three quino- 
line rings. Note that we have broadened the DOS peaks, 
leading to a shift to lower energies. The LUMO consists 
only of atomic p-type orbitals, no other type is present. 
One can observe that electron doping commences with 
initially filling the empty p z levels of N3. Hence, the mo- 
ment on N3 increases at first steeply with doping. Due to 
the small asymmetry and the structural relaxation, the 
empty p x and p y states of Ni and N2 start to be filled 
more slowly. Note that these are equivalent to the N3 
p z orbitals by rotation. Electron doping of more than 
0.5 e leads to an increase of the moments on all three ni- 
trogens at the same rate. For the narrow p z orbitals 
spin-polarized filling is more favorable than equal fill- 
ing of both spin-orbitals (Pauli principle) . A somewhat 
larger electron filling would lead to an equal filling of both 
spin-orbitals on one nitrogen and thus to a non-magnetic 
state, but this does not happen, because the unoccupied 
p z orbitals of the other two nitrogens are preferably be- 
ing filled first, which causes a small moment on Ni and 
N2 as well. A ferromagnetic coupling between the three 
spin moments stabilizes finally a net spin-polarization on 
the whole molecule. 

To conclude, our DFT calculations reveal that electron 
injection in Alq3 leads to spontaneous spin-polarization. 
The predicted spin-polarization reaches to 1 hb per in- 
jected electron per molecule. The spin-polarization is 
found to originate from the preferential spin-polarized 
filling of empty N p z orbitals in the LUMO, in com- 
bination with a parallel coupling between the spins on 
the quinolincs. This discovered induced spin-polarization 
of the normally non-magnetic Alq3 molecule can be ex- 
pected to be critically related to the extremely long spin 
relaxation times observed in Alq 3 Q as well as to be 



essential for the observed large magnetoresistance [l[. 
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